Why the two-pulse photon echo is not a good quantum memory protocol 
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We consider in this paper a two-pulse photon echo sequence as a potential quantum light storage protocol. It 
is widely believed that a two-pulse scheme should lead to very low efficiency and is then not relevant for this 
specific application. We show experimentally by using a Tm''"'" : YAG crystal that such a protocol is on contrary 
very efficient and even too efficient to be considered as a good quantum storage protocol. Our experimental 
work allows us to point out on one side the real limitations of this scheme and on the other side its benefits 
which can be a source of inspiration to conceive more promising procedures with rare-earth ion doped crystals. 
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PACS numbers: 03.67.-a, 32.80.Qk, 42.50.Md, 42.50.p, 42.50.Gy, 42.50.Md 



I. INTRODUCTION 

The prospect of quantum light storage in solids motivates 
us to reconsider the interaction of light and matter at the sin- 
gle quantum level. Historically coherent transient phenomena 
appeared very active because they are primarily based on two- 
level system absorption. The community naturally focused 
on rare-earth ion doped crystals (REIC) instead of atomic va- 
pors because the storage time can be long as well and is not 
limited by the atomic diffusion whatsoever Recent propos- 
als are of course inspired by previous realizations of classical 
light memories [H^ and all-optical processing more generally 
02]. In this lineage recent progresses toward quantum stor- 
age are involving common physical ingredients. Because of 
the large inhomogeneous broadening, any light retrieval is in- 
timately related to a dipole rephasing. Surprisingly the very 
well known conventional photon echo has not been considered 
as a promising quantum protocol despite impressive realiza- 
tions in the classical domain. A major identified drawback is 
the low efficiency of the process ylHl- This has been a com- 
monly admitted idea since the pioneering studies of the 60's 
llH] . Although we roughly adhere to this statement, we would 
like to point out that high retrieval efficiency was already ob- 
served in a specific regime of two-pulse photon echo (2PE) 
10, 1?]], which might pave the way to quantum light storage. 
Our experimental work consists of a clear observation of the 
large predicted efficiency in this regime we first need to pre- 
cisely defined. Our analysis is also stressing clearly what are 
the physical ingredients which lead to this result. Our study is 
definitely placed in the prospect of a quantum memory. With 
respect to other protocols, we finally clarify the advantages 
and drawbacks of this technique that should be considered as 
a general tool for coherent manipulations. 

As mentioned before, the investigation of classical light 
storage largely paved the way toward their quantum equiva- 
lent. REIC have shown interesting processing capabilities es- 
pecially with all-optical control |2|, [8|]. As derived from the 
conventional photon echo, these techniques are based on an 
optical manipulation of the coherences. Experiments largely 
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benefit froin the agility of the laser controlling the crystalline 
processor [9l]. This convenience would be still appreciable 
for manipulation at the quantum level. The 2PE time-to- 
bandwidth product properties should also be emphasized. In 
the 2PE process, this parameter, critical for information pro- 
cessing applications, is not limited by the memory opacity, 
in contrast with the most promising quantum storage proto- 
cols involving REIC, namely the "stopped-light" approach 
ifioi ITTI] or the "controlled reversible inhomogeneous broad- 
ening" (CRIB) procedure OHlllIllTl. Finally, contrary 
to the above storage protocols, doesn't require any initial state 
preparation. A spectral selection within the inhomogeneous 
broadening is in a sense build-in because of the selective ex- 
citation of the first incoming pulse that we define as the sig- 
nal. The 2PE has the singular advantage to rephase a ran- 
dom distribution of level shifts without any assumption on the 
source of inhomogeneity. By clarifying the characteristics of 
the 2PE, we aspire to a deeper understanding of the atomic 
coherences optical manipulation. 



Since the observation fs*] and the interpretation fW\ of the 
photon echo, the quantitative comparison with the observed 
efficiency has been widely studied. At the basis of data pro- 
cessing apphcation, the interest has been renewed relatively 
recently Oil 11211 ■ A realistic approach of the problem usually 
requires a numerical resolution of the Bloch-Maxwell equa- 
tions jTsl. This gives a solid interpretation of the experimen- 
tal data [171. As a matter of fact investigating the 2PE for 
quantum light storage is strongly simplifying the problem. As 
compared to the canonical 2PE where a 7r/2 excitation pulse is 
followed by a vr rephasing pulse, memory-like version of this 
scheme would first involve a very weak signal pulse. In that 
very specific case as pointed out by Tsang et al. |,19,1, one can 
derive an analytic solution for the efficiency. In this paper, 
we will first deduce these equations from a simple physical 
interpretation of the scheme. We will clearly specify the un- 
derlying assumptions to be verified in practice. We will then 
compare these calculations to the experiment. A detailed anal- 
ysis of our observations allows us to conclude and place this 
work in the context of a quantum memory by comparison to 
other storage protocols. 
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II. EFFICIENCY OF THE PROTOCOL 



A. Calculation in tlie weak signal limit 



This subject has been covered by a wide range of literature 
10, S [13, [IE [isll- Nevertheless it is relatively easy to derive 
these equations based on simple physical arguments. The in- 
teraction of light pulses with our medium is well described 
by the Maxwell-Bloch equations assuming the slowly varying 
amplitude and the rotating wave approximations. 



dtn{ujab]z,t) = -Av{uJab;z,t) 

dtv{uJab; Z, t) -^{z, t)w{uJab; Z, t) + l^n{uJab\ Z, t) 

dtw{uJab;Z,t) ^fl{z,t)v{uJab;Z,t) 

(1) 

where A = ujab — is the detuning, f2 is the Rabi frequency 
of the field under consideration and (u, v, w) the three com- 
ponents of the Bloch vector The decay of the coherences and 
the population is assumed to be negligible. We have dropped 
the usual term dt^{z,t) because realistically the spatial ex- 
tension of the pulse is always much longer the length of our 
crystal. 

To describe the broadest range of situation a numerical res- 
olution of the system is usually necessary |18]. This is not 
our approach. A sketch of the time sequence is depicted in 
FigdJand looks like any 2PE sequence. Nevertheless within 
the prospect of quantum storage, the signal and the echo are 
assumed to be weak. This greatly simplifies the description 
[|20|] essentially because these two fields do not modify the 
population difference w{LUab] z, t) that is not time dependent 
anymore. This is the small area approximation where the 
Maxwell-Bloch system can be linearized 1I20I1 . 
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FIG. 1: Outline of the time sequence. The signal is first absorbed 
in the medium. After a time ti2 a rephasing pulse of large area A 
induces a build-up of the coherence at the time 2ii2 and give rise to 
the photon echo. We address the problem in the weak signal and the 
echo regime, which shall be satisfied in quantum field conditions. 



If the three pulses are well separated in time, one can con- 
sider them independently. The signal, the rephasing pulse and 
the echo are respectively centered on t = 0, ti2 and 2ti2. The 
incoming signal and outgoing echo Rabi frequencies are re- 
spectively denoted S and £. The propagation of S is simply 
described by an absorption law if the atoms are initially in the 
ground state w{ujab'-, z, —00) = w{uJab'-, z, 0) = —1 |20]: 

dzS{z,t) = w{wab;z,Q)^S{z,t) = ~^S{z,t) (2) 

The weak echo is expected to behave in a similar way, ex- 
cept the medium has been previously excited and modified by 
the signal and the rehasing pulse. Therefore the echo equation 
reads as: 



dz£{z,t)^ w{uJab;z,2ti2)j£{z,t) 

~2tT '^ab'^-E{0Jab;Z,t) 



(3) 



The coherence vg, resulting from interaction with the first two 
pulses, evolves freely within the time interval 2ti2. The 

population w{LUab', z, 2^12) has been affected by the rephasing 
pulse. We assume the rephasing pulse is much shorter than S. 
Therefore the population is uniformly modified by the rephas- 
ing pulse all over the spectral interval initially excited by S. 
We will see in|lII]how this constraint is treated experimentally. 
A second benefit of this assumption is that the rephasing pulse 
can be consider as instantaneous (time ^12) and then uniformly 
modifying the coherences. The physical interpretation of Eq|3] 
is based on the generation of £ by the macroscopic dipole that 
builds up when the coherences ve get phased together. The 
growing field propagates through the medium characterized 
by the uniform population difference w(aja6; z, 2^12). 

We aim at reducing Eq[3]by expressing the atomic quanti- 
ties in terms of the optical fields only. We shall be left with an 
equation of propagation for £. The echo efficiency will be de- 
duced from the solution of this equation. We first have to track 
the excitation of the coherences by the signal, then their mod- 
ification by the rephasing pulse and finally their free evolution 
toward the echo emission. The problem is addressed locally, 
at position z. The signal S excites the atoms that initially all 
sit in the ground state: 



u{u;ab]Z,t) +iY{uJab;Z,t) 



i exp (iAt) / S{z, t) exp (— i At) dr 

J —00 



(4) 



At a certain time t between and ti2 the signal field is off. 
We then recognize the Fourier transform of S written S, the 
exp (zAi) accounts for the free evolution during this interval. 

\i{ujab; z, t) + iv(cjQ6; z,t) ^ i exp (iAt) S {z, A) (5) 

This expression represents the evolution of the coherence after 
the initial absorption process and before the rephasing pulse. 
Next, we calculate this strong pulse effect on the coherences 
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(how ve is related to v) and on the population, which modifies 
the echo propagation (Eql3]). This can be done analytically 
by integrating the Bloch-Maxwell equations (see ref. lO). 
Nevertheless the results are relatively intuitive at the end and 
can be derived from simple physical ingredients. This can be 
done first of all in the specific case of a 7r-rephasing pulse, 
the more general case of an area A{z) for the strong pulse 
can be solved by introducing "by hand" geometrical factors. 
All over the spectral interval excited by the signal pulse, the 
rephasing pulse is assumed to behave as a vr-pulse. This cor- 
responds to a brief pulse assumption. The experimental ful- 
fillment of this condition will be addressed in |III] A vr-pulse 
simply drives the Bloch vector by a rotation of tt around an 
equatorial axis. On the one hand, along the population axis it 
corresponds to an inversion from —1 to w{u!ab', z, 2^12) — 1 
at the time ti2. On the other hand, it transforms the coher- 
ences v(wq{,; z,ti2) -v{ujab; z,ti2), u(wa6;z,ii2) stays 
the same (complex conjugation of u + iv). Right after the 
rephasing pulse, Eq. |5]becomes 

UE(Wab; Z, t) +iVE(t^afc; 2, t) = 

i exp {iA {t - ti2)) S* (z, A) exp (-iAti2) 

(6) 

where the complex conjugation sign * accounts for the rephas- 
ing transformation. The coherences are freely evolving af- 
ter ti2- One can now write the propagation equation of the 
echo (Eq|3]i by including the direct influence of the signal on 
Ve i^ab ; 2 , i ) . The signal field being assumed to be a real num- 
ber, one recognizes the time-reversed signal S{z, 2ti2 — t) 
whose z dependency is given by an absorption law 2ti2 — 
t) = 5(0, 2ii2 - t) exp (-Q!z)(Eq|2li 



dz£{z,t) 



+ -£{z,t) 



aS{z,2ti2-t) 



(7) 



The signal field acts as a source and generates the echo that 
propagates in an inverted medium. This gives the equation 
of propagation for a 7r-rephasing pulse. It is now rather 
easy to account for an imperfect rephasing. More gener- 
ally, a yl(z)-area strong pulse drives an j4(z)-rotation of the 
Bloch vector The population is not fully inverted anymore: 
w{ujab', z, 2ti2) — — COS A (z). The rotation of the coherences 
is also incomplete and limited to [1 — cos {A (z))] /2 of its 
maximum value. This factors are purely geometrical and are 
interpreted as projections on the Bloch sphere. We finally get 
the general analytic expression for the efficiency. This expres- 
sion has been previously derived by Tsang et al. fl^, Eq.(40)] 
by integrating the Bloch-Maxwell equations. Here we simply 
focus on the two crucial stages, the absorption of the signal on 
one side and the re-emission of the echo on the other side. The 
rephasing pulse in between is interpreted as an instantaneous 
manipulation of the Bloch vector. 



d^£{z,t)= -cos(A(z)) f£(z,i) 
1 - cos {A (z)) 



(8) 



aS{z^ 2ti2 — t) 



There is an underlying assumption here: the rephasing pulse 
is very brief and is then fully covering the spectral range of ex- 



citation. As a consequence the echo is not deformed as com- 
pared to the signal, it is only time reversed. The ^(z)-rotation 
on the Bloch sphere is also uniform and doesn't depend on 
A. The z-dependency of the area accounts for the propaga- 
tion of the strong pulse itself. This is usually a complicated 
problem but in that case we are only interested in the propaga- 
tion of the area. The result is remarkably simple and is given 
by the Area Theorem of McCall & Hahn |21]. We don't have 
to know the exact temporal shape through the propagation be- 
cause the area is the relevant quantity for the rephasing pulse 
and is simply given by 



d,A{z) = -- smA{z) 



(9) 



It can be solve analytically for a given A{0). A straightfor- 
ward integration of Eq|8] allows us to calculate the retrieval 
efficiency 77 as a function of the optical thickness aL where L 
is the length of the medium: 



rj{A{0),aL) 



£{L,t) 



5(0, 2^12 -t), 

2 sinh {aL) 
1 + exp [aL] cot2 {A (0) /2) 



(10) 



For different optical thickness aL, we plot the efficiency as a 
function of the rephasing pulse area A (0) in Fig|2] We ob- 
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FIG. 2: Efficiency of the 2PE as a function of the rephasing pulse 
area A (0) for different optical thickness otL. Inset: for larger otL, 
the efficiency is much larger than unity at the maximum rephasing. 

serve in that case that the efficiency is strongly depending on 
the optical thickness. When it is low, the efficiency is weak es- 
sentially because the signal is poorly absorbed. The efficiency 
is then a sinus-like function and the rephasing area directly 
accounts for an imperfect rotation on the Bloch sphere. On 
the other side, at large optical thickness, the efficiency can be 
much larger than unity but only on a narrow window around 
a TT. This specific situation is certainly the most interesting 
because the signal in fully mapped into the medium (large op- 
tical thickness) and the efficiency is exceptionally large. The 
TT -rephasing pulse have very particular properties that can be 
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interpreted independently to physically understand this behav- 
ior. 



B. Specificity of tlie 7r-repliasing pulse 

The case of an exact 7r-rephasing pulse is relatively straight- 
forward. More fundamentally it allows us to understand the 
expected efficiency at large optical thickness. We simply de- 
rived the propagation equation O for the echo by assuming 
that locally, at position z, the area is exactly tt. Now we have 
to examine the 7r-pulse propagation inside the medium. This 
is a very specific situation. According to Eq|9]a 7r-pulse pre- 
serves its initial area throughout the medium. Even if the en- 
ergy is absorbed, the area is conserved as the pulse stretches 
temporally (tS]. Roughly speaking, if the pulse is elongated 
by a factor r, the amplitude (Rabi frequency) is reduced by ?■ 
to conserve the area and the energy decreases by a factor r. 
This alteration is a pure coherent propagation effect. 

The propagation equation [7] is then valid at any position z 
and easily gives the efficiency. This expression is consistent 
with our general formula Eq. [TO) 

?7 (tt, aL) = [exp (ai/2) - exp {-aL/2)f (1 1) 

At large optical thickness, the efficiency is much larger than 
unity and grows exponentially. This is relatively counter- 
intuitive. The echo efficiency is generally observed to be low, 
which is usually assigned to absorption. As mentioned before 
the TT-rephasing pulse retains its area along the propagation. 
In other word, the medium is completely inverted: the echo 
is emitted in an amplifying medium. This explains why the 
echo is gaining exponentially. Practically, the assumptions we 
made, such as 1 -dimension, infinite plane wave geometry, will 
be difficult to satisfy. Any divergence from the ideal theoreti- 
cal frame shall affect the echo efficiency. 

The TT-pulse propagation is not only unusual, it is also a 
singular solution of the Area Theorem. The vr-solution is 
indeed not stable because any area slightly lower (or larger) 
than TT will decrease (or increase resp.) toward (or 27r resp.) 
ll2ll l22h. Even so, the pulses with an area close to tt can prop- 
agate deeper inside the medium than a weak pulse. To see that 
we plot in Fig. [3]the penetration depth Lp at which the incom- 
ing area is divided by 1/ -/e (the curve is symmetrized around 
TT to account for the deviation toward 2tt of pulses larger than 
t:). 

For small areas, this length is a^^ as expected ll20ll . As 
one gets closer to tt, the pulse can propagate deeper inside the 
medium. For example, if the incoming area is controlled at 
the 1 % level (inset Fig. O, the penetration depth is larger than 
7a~^. The curve is very narrow about tt which shows a high 
sensitivity for the propagation. Realistically a well-controlled 
TT-pulse should induce a population inversion much deeper in- 
side the medium than the absorption length a^^. If this depth 
is larger than the optical thickness of the medium, it will be 
fully inverted. However, the deeper the 7r-pulse propagates 
through the absorbing medium, the more it is stretched, since 
it has to keep a constant area while losing energy. The pulse 
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FIG. 3: Penetration depth Lp of a strong pulse in unit of as a 
function of the incoming area. Inset: Penetration of pulses with area 
of TT plus or minus 1%. For small area pulses, it is simply given by 
an absorption law Lp — . 



bandwidth shrinks accordingly, making the pulse act as a tt- 
pulse on a reduced spectral interval |22]. 

Based on this analysis, we expect to observe two remark- 
able qualitative features. The efficiency should be very high 
at large optical thickness and strongly depending on the area 
of the rephasing pulse. 



III. EXPERIMENTS 

A 2PE experiment can be performed in any system were a 
transient phenomena can be observed. Here we use a thulium- 
doped yttrium aluminum garnet (YAG), cooled down to cryo- 
genic temperature. The long optical coherence time makes 
it particularly attractive for quantum storage application. We 
will now briefly describe the experimental set-up and focus on 
the precautions we take to satisfy the assumptions introduced 
previously. 



A. Experimental set-up 

Our 0.5% Tm'^+: YAG crystal is immersed in liquid helium 
at 1.4K. The coherence time of the ^H6(0) to ^H4(0) transi- 
tion is typically T2 ~ in these conditions. The crys- 
tal is oriented and cut in order to propagate along the [lIO] 
direction. Along this axis, the length is 5mm and the opti- 
cal thickness aL = 5. The laser polarisation is parallel to 
[111] to maximize the Rabi frequency ||23|] . The laser system 
is operating at 793 nm, stabilized on a high-finesse Fabry- 
Perot cavity with the Pound-Drever-Hall technique (200Hz 
over 10ms) [i24j]. The laser is split in two independent beams. 
Temporal shaping is achieved by two acousto-optic modula- 
tors (AOM) controlled by a dual channel arbitrary waveform 
generator (Tektronix AWG520). Both beams are injected into 
two single mode fibers. Before recombination on a beamsplit- 
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ter, we use expanders to independently manage their waists 
inside the crystal. After the sample, the signal is collected in 
a single mode fiber terminated by a photodiode. 

The signal is supposed to mimic a weak quantum field, so 
this pulse should verify the small area approximation. The sig- 
nal beam is in practice much weaker than the rephasing one. 
At the maximum, there is a 36dB power difference between 
the two. More precisely, the signal area is kept constant at 9% 
of TT and we vary the rephasing pulse area A (0) from to Stt. 

Keeping the rephasing pulse significantly briefer than the 
signal pulse is the most stringent condition we have to satisfy. 
This is required to maintain a uniform coverage of the rephas- 
ing process over the signal excitation bandwidth. To do so, 
we use gaussian-shaped signal (duration 2.1/is) and rephasing 
pulses. On the one hand a gaussian pulse is spectrally nar- 
rower than a rectangular pulse with the same duration. On 
the other hand we observed a gaussian 7r-pulse undergoes 
less temporal streching than a gaussian one after propagation 
through an absorbing medium. We set the rephasing pulse 
K 2.5 times shorter than the signal. This value is slightly fluc- 
tuating depending on the rephasing amplitude value. Chang- 
ing the AOM driving power marginally impacts on the pulse 
shape. A much shorter pulse would be preferable but we are 
limited by the available power (few milliwatts) to ensure a 
significant area in a reasonable time. 

One last point we ignored so far is the transverse dimen- 
sion of the beams. To be consistent with the 1 -dimension 
theory, the power of the rephasing beam should be constant 
over the spatial extension of the signal. The signal beam waist 
(17/im) is then chosen to be 2 times smaller than the waist of 
the rephasing beam (35/^m). This is the same overlap argu- 
ment we used in the spectral domain. 



B. Results 

We perform a 2PE experiment in the beforementioned con- 
ditions (see Fig. |4]i. We pay special attention to an accurate 
calibration of the rephasing pulse area. We indeed first per- 
form an optical nutation experiment to evaluate the exact Rabi 
frequency of the pulse. Comparing the signal and rephasing 
beam intensities, we estimate the area of the signal which is 
confirmed to be weak (O.OOvr). 

As expected, we observe an echo at delay ti2 after the 
rephasing pulse (see Fig. |4]). We also carefully calibrate the 
efficiency. To do so, we shift the laser far from the absorption 
line (a few cm^^). The corresponding intensity level repre- 
sents the 100% reference line on the measurement detector. 
According to the efficiency definition given by Eq. [TO] the 
echo is assumed to exhibit the time-reversed temporal shape 
of the signal. This is not exactly the case experimentally as we 
shall discuss later So we define the efficiency by comparing 
the maxima of the two pulses: 



x10" 



ryexp (A(0)) = 



/ maxt (g(L,t)) 
\maxt (5(0, t)) 



(12) 



By varying the incoming rephasing area, we obtain 77cxp as a 
function of A (0) (see Fig. |5]l. 
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FIG. 4: (color online) Time sequence for Llyr-rephasing pulse. The 
signal and the echo being very weak, we magnify their scale by a 
factor 150 (solid red line). The intensity is given as the square of the 
Rabi frequency. The signal and the rephasing pulse have a duration 
of and respectively (rms width of the gaussian). We also 
represent the rephasing pulse after transmission (dashed line in Arb. 
units). 
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FIG. 5: Efficiency of the echo retrieval as a function of the rephasing 
pulse area A (0). The maximum (circle) is obtained for A (0) = 
I.Itt where the efficiency is 92%, this corresponds to the situation 
depicted in Fig. |4] 



The main source of uncertainty is due to the alignement and 
the spatial overlap of the beams in the crystal. To quantify it 
and derive error bars, we estimate the typical intensity vari- 
ation of the rephasing beam (waist 35/im) over a length cor- 
responding to the signal waist (17/im). This is simply given 
by the direct comparison of two gaussian curves (11% in that 
case, which gives the horizontal error bars in Fig. |5]l. 
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C. Discussion 



The experimental results are qualitatively in good agree- 
ment with the expected efficiencies in Fig. ID The main fea- 
tures are indeed observed. First of all the maximum efficiency 
is obtained for a l.l7r-area, which is consistent with tt, within 
the error bar. We also observe a reincrease of the efficiency 
close to Stt after a minimum at 2-k. We certainly predict a 
maximum at any odd number of tt, but we are experimentally 
limited by the available laser power. Secondly, the curve is 
peaked in the vicinity of vr and cannot be fitted by a sinus- 
shape oscillation. This is also expected (see Fig. |2]i and is 
due to the large optical thickness of the sample aL ~ 5. 
Finally, the maximum efficiency is relatively high 0.92 for 
A{Q) — 1.1 TT (circle in Fig. |4]l. Although far below the 
rj (tt) = 146 predicted value (Eq. fTTT i. this result demonstrates 
a highly efficient 2PE. 

We can invoke many reasons to explain the discrepancy be- 
tween the measured and the predicted efficiency values, (i) 
The first obvious one is the total duration of the time sequence. 
The echo is indeed decaying exponentially because of the co- 
herence lifetime T2 — 50fis which has been completely ne- 
glected in our treatment. In our case ti2 = 25/iS so the echo 
is reduced by a factor cxp(4ti2/T'2) = 7.4. Without this de- 
cay, the efficiency would be much larger than unity. In section 
llVlwe shall see why ti2 is chosen to be long for this experi- 
ment, (ii) Another limitation is certainly due to the duration 
of the rephasing pulse. With this pulse 2.5 times shorter than 
the signal, the incoming pulse spectral overlap is rather good. 
However, propagation through the sample strongly stretches 
the rephasing pulse, as expected from discussion in Sec. IIIBI 
and as observed in Fig. |4] The pulse cannot be considered 
as much briefer than the incoming signal all the way through 
the sample. An observable proof of this effect is the retrieval 
time of the echo (Fig. |4|i. The retrieval should be centered on 
t — 2ii2 — 50/is. We clearly see that the echo is delayed by 
few microseconds. Indeed, because of stretching the rephas- 
ing pulse is no longer centered at a delay ti2 from the signal. 
The retrieval time is shifted accordingly. Since one of the as- 
sumption of our model is not fully verified, we then expect an 
efficiency reduction, (iii) In the spatial domain, the same argu- 
ment is also valid. The signal is tightly focused (17/im) to en- 
sure that its waist is smaller that the rephasing beam. The as- 
sociated confocal parameter is typically two time shorter than 
the crystal length. So the rephasing beam does not overlap the 
signal uniformly all along the propagation. In other words, the 
rephasing area is not constant in the transverse direction. This 
should reduce the efficiency and broaden the peak around tt 
(convolution effect). 

Based on this analysis, we believe our model contains all 
the physical ingredients to explain qualitatively the experi- 
mental results. We have given three probable explanations to 
interpret the quantitative discrepancy with the predicted val- 
ues. 



IV. RELEVANCE FOR QUANTUM MEMORY 
APPLICATION 

Our experiment has been performed in the classical domain 
using weak small-area pulses. It tells us however what should 
be the limitations in the quantum domain. 

An obvious one is already present in our experiment. Since 
the rephasing pulse stretches while propagating through the 
sample, it gains a trailing tail that is not negligible as com- 
pared to the echo amplitude. As we can see in Fig. |4] at 
t = 45/iS the pulse tail is falling slightly before the echo 
comes out. That's the reason why we cannot make ti2 shorter, 
otherwise the echo would be submerged. With only few pho- 
tons in the signal, this effect would be disastrous. As already 
discussed, the strong pulse distortion is not an artefact. This is 
a coherent propagation effect |22], thus a fundamental limi- 
tation. This should not be confused with the noise induced by 
the fluorescence, which will be another limitation at the few 
photons level. 

The rephasing process is inherently associated with a pop- 
ulation inversion. The decay has been neglected in our model. 
In practice the medium excitation will be followed by sponta- 
neous emission. The fidelity of the 2PE as a quantum mem- 
ory protocol is fundamentally limited by fluorescence. This 
can akeady be understood within the framework of the Dicke 
model i25ll . i.e. without taking propagation effects into ac- 
count. Consider an ensemble of two-level systems, where the 
two states are denoted \g)k and \e)k for the fc-th system. Note 
that we are interested in the case where the transition energy 
for the g to e transition is slightly different for different sys- 
tems (inhomogeneous broadening). 

We will compare the case where the input to be stored is a 
single photon to the case where there is no input (i.e. where 
the input state is the vacuum). The initial state of the atomic 
ensemble is 1-0") = \g)i\g)2---\g)N- For a vacuum input, this 
state remains of course unchanged. It is then transformed to 
lip^) — |e)i|e)2...|e)Ar by the vr-pulse. For a single photon 
input, absorption of the photon creates a state of the form 

= ^ (|e)i|g)2...|5>A^ + - + |.g>i|g>2...|e)jv) (B) 
V iV 

which contains a single atomic excitation. (Propagation ef- 
fects would lead to the coefficients of the N terms not being 
all the same.) The 7r-pulse transforms this state into 

IV-^-') = ^ (|<?)i|e)2...|e)^ + ... + |e)i|e)2...|<?)^) 
V iV 

(14) 

which has iV — 1 atomic excitations. 

In an inhomogeneously broadened system, the various 
terms in Eqs. ( I13I14I ) will acquire different phases depend- 
ing on the transition energies of the various atoms. However, 
at the time of the echo all terms will be in phase. Emission 
from the state gives rise to the echo signal corre- 

sponding to a single-photon input, whereas emission from the 
state IV-'^) corresponds to a vacuum input and thus defines 
the noise background due to fluorescence. The photon emis- 
sion probability for a state \^) is proportional to ||J_|^)|p, 
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where J_ = X]/c=i l5)fe(sU- This is due to the fact that the 
interaction Hamiltonian between the atomic ensemble and the 
relevant mode a of the electro-magnetic field (corresponding 
to emission in the direction of phase matching) is proportional 
to J- + h.c. 

Following Ref. iQ it is easy to see that 1 1 J_ I ? = 

2{N — 1) and || J_|i/i^)|p = TV. As a consequence, the prob- 
ability to emit a photon at the echo time is only twice as large 
for a single-photon input as for no input at all, correspond- 
ing to a signal-to-noise ratio of one. This severely limits the 
achievable fidelity of quantum state storage. 

Finally, from what we have shown is this article, we can 
conclude that the efficiency is actually too high for quan- 
tum memory application. We indeed observed a maximum 
92% retrieval but the efficiency can be much larger than unity 
with optimized conditions [6J. This amplification due to the 
medium inversion is precisely a propagation effect that is not 
considered in the beforementioned Dicke model. This is a 
key ingredient to interpret our experimental results. In quan- 
tum optics terms, the statistics of the field will be modified: 
for one photon coming in, more than one will come out. Be- 
cause the medium is inverted, it actes as a gain medium and 
modifies the quantum field and then again reduces the fidelity 
ll26ll . In that sense, the 2PE is also too efficient to be a good 
quantum protocol. 

Those reasons are three fundamental limitations that we ex- 
pect in the quantum regime. Even if the 2PE suffers from 
drastic drawbacks, it should be considered with attention. It 
is not only a historical example that helps us to understand 
rephasing phenomena. It has the unique ability to rephase 
atoms with randomly distributed level shifts, whatever the dis- 
tribution structure Iil41 . This a major difference as compared 
to CRIB (Controlled Reversible Inhomogeneous Broadening) 
1(121 [T3I [3, [Htl where the first initial step is a spectral tailor- 
ing of the distribution. There is not preparation of the system 
in the 2PE. The feat performed by the CRIB protocol is its ca- 
pability to produce a rephasing of the coherence without any 
population inversion. As a consequence the equations of the 
2PE (Eq. |7]) and the CRIB lf27ll are remarkably similar except 
that a minus sign accounts for the population inversion in the 
propagation equation. One can finally wonder if an optical 



manipulation would achieve a rephasing in the ground state as 
the CRIB does. This is a priori not possible because even a 
complex optical sequence will be decomposed with rotations 
on the Bloch sphere, on contrary the CRIB protocol can be in- 
terpreted as a planar symmetry (detuning sign reversal). These 
two are then intrinsically and fundamentally different. 

We have here listed the limitations of the 2PE when consid- 
ered as a potential quantum storage protocol. We pointed out 
the pulse deformation than can be a technical issue when us- 
ing strong pulse. The two other limitations are directly and 
fundamentally related to the medium inversion induced by 
the optical rephasing operation. On one side the spontaneous 
emission will produce a noise comparable to the retrieved sig- 
nal and then deteriorate the storage fidelity. On the other side 
the inversion will make the medium amplifying, which mainly 
explains the large predicted and observed efficiency. A larger 
than one efficiency is also associated with a fidelity reduction 
for quantum fields. 



V. CONCLUSION 

In this article, we study the 2PE efficiency within the con- 
text of quantum storage. In this framework, we experimen- 
tally observe large efficiencies that are well explained by a 
simple model. Our calculations are based on a physical analy- 
sis of this specific situation. The experimental set-up has been 
devised to verify the underlying assumptions of the model. In 
2PE, rephasing goes along with population inversion. This 
is a crucial ingredient of this protocol. The emitted echo is 
then widely amplified and can be stronger than the incoming 
signal. We have observed this effect. We finally conclude 
by analysing the potential extension of this work at low light 
level. We have pointed out the inherent limitations of the pro- 
cess. 

By clarifying the physics involved in the very well-known 
two-pulse photon echo, we more generally tackle the problem 
of using strong light pulses for rephasing purposes. Our study 
should then be considered as a tool for the conception of new 
quantum storage protocols. 
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